Analytic behavior of the QED polarizability function at finite temperature 
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We revisit the analytical properties of the static quasi-photon polarizability function for an elec- 
tron gas at finite temperature, in connection with the existence of Friedel oscillations in the potential 
created by an impurity. In contrast with the zero temperature case, where the polarizability is an 
analytical function, except for the two branch cuts which are responsible for Friedel oscillations, 
P5 ' at finite temperature the corresponding function is not analytical, in spite of becoming continuous 

' everywhere on the complex plane. This effect produces, as a result, the survival of the oscillatory 

behavior of the potential. We calculate the potential at large distances, and relate the calculation 

■ to the non-analytical properties of the polarizability. 

<^ . 

I. INTRODUCTION 

(-H I The potential created by a static ionic impurity in metallic alloys has been considered by many authors (see, for 
example [1, 2]). At large distances, the potential shows an oscillatory behavior, damped as negative powers of the 
distance r, a phenomenon which is known as Friedel oscillations At zero temperature, this kind of behavior has 
. been associated to the existence of the Kohn singularity Q in the quasi-photon polarizability, induced by the sharp 
edge of the degenerate electron distribution. At the threshold of the interaction, i.e. at the threshold of electron-hole 
, creation, the momentum of the quasi-photon is equal to 2p/, the diameter of the Fermi sphere, with pj the Fermi 
^ momentum of the electrons. It must be noticed that the quasi-photon self-energy presents a singularity at this point. 
C\| , The presence of a singularity allows to obtain, with the help of Lighthill's method, the asymptotic form (r oo) of 
1 the potential, as an expansion in terms of the form cos{2p fr) and sin(2p/r), damped as negative powers of r, and 

■ enhanced by powers of log(4p/r) f^. 
Alternatively, one can perform an analytical continuation of the quasi-photon polarizability to complex values of 

the transferred momentum q, and do the integration by deforming the circuit in the complex plane. Aside from the 
"Debye pole" of the quasi-photon propagator on the imaginary axis, giving an exponentially damped contribution, 
the polarizability has two branch cuts starting at q — ±2pj, which are responsible for the long-distance oscillatory 
behavior [l|, Q . In other type of plasmas, one can also find additional complex poles in the propagator, giving raise 
to an exponentially damped oscillatory contribution frj. 

Let us now consider the situation at finite temperature T. Since the electron distribution is spread out, from 
energy and momentum conservation in the collisions between the electrons and soft quasi-photons we expect that 

■ screening becomes more effective than in the zero-temperature case. This is, in fact, the case and one obtains that 
the oscillations are still present, although they are damped as temperature increases Of course, this is a desirable 
result, showing that the limit T = is not pathological. 

We would like to understand this aspect from a mathematical point of view. On one side, the above discussion, 
based on the Kohn singularity, does not hold, since we know that the singularity disappears at T ^ On the 
other side, one can be tempted to extend the integration to the complex q -plane. If the T > polarizability is 
still an analytical function, one has to face these possibilities: either the function has discontinuities on the complex 
plane (i.e. branch cuts, as in the T = case) or poles , or both. As we discuss later, the first possibility does 
not appear: the polarizability becomes a continuous function at non-zero temperature. On the other hand, the 
appearance of additional poles, at an arbitrary small temperature, which were not present at T = 0, would indicate 
a strong anomaly, as branch cuts should be suddenly transformed into complex poles, a possibility which looks too 
exotic and would point towards a singular behavior with temperature. As discussed above, this is not the case. Then, 
how can we account for the oscillations? 

As we show in this paper, the explanation lies on the fact that, at finite temperature, the polarizability is non- 
analytical over the whole complex plane. In fact, this function can be understood as a superposition of a family of 
functions, each one having discontinuities at g = ±2p. Each one of the functions on this family gives an oscillatory 
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result and the sum, for not too large temperatures, is still oscillatory. As temperature increases, the range of values of 
p that effectively contribute is enlarged, leading to a destructive interference. For this reason, oscillations are damped 
with temperature. 

This paper is organized as follows. The analytical properties of the polarizability function, both at T = and 
T > 0, are discussed in sections 2 and 3. In section 4 an expression for the potential is obtained. Section 5 contains 
a brief discussion and concluding remarks. 



II. THE POLARIZABILITY IN THE COMPLEX PLANE 

The potential created by a static ionic impurity in an electron gas can be calculated from: 



V{r) = —I{r) 

I ir) - r dq ^ r ii^, (1) 

J-oo nq) J-oo nq) 



where ll!: 



%)-g'-Ax^('Z), (2) 

A = 47re^, e is the electron charge, and x^i^) is the polarizability at temperature T within the RPA (random phase 
approximation), which can be written as: 

/"OO 

X^{q) = / dpgip)x"{q,p). (3) 



Here, 

X°{q,p) = m ' J"^^"""^' (4) 

is the T = polarizability. m is the electron mass, and 
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9iP)--^ (5) 



where / is the Fermi-Dirac distribution for the electrons 



l + exp(^^ 

We consider a non-relativistic plasma. Therefore, in the latter formula e {p) — ^ . Finally, is the chemical 
potential. Under these conditions, we have: 



As mentioned above, Eq. ^ corresponds to the RPA level. A further improvement of this approximation can be 
introduced via a local field correction. For zero temperature, simple analytical modelizations are given e.g. in 
However, to our knowledge no such analytical models exist at finite T. For this reason, since we only intend to give 
a qualitative explanation for oscillations at T > 0, we stay at the RPA level. Of course, in order to obtain more 
accurate results, one would need to go beyond this approximation. 

The usual way to calculate the integral in Eq. ([1]) is by performing an analytical continuation to the upper complex 
half-plane: 

m^Imfdq'-^^, (8) 
Jc v{q) 

where C is now a circuit contained on the upper half-plane. In doing so, we need to investigate the analytical properties 
of ^^{q) (see Eq. [5]) in the interior of C. These properties are defined by Eq. ([3]). 
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III. ANALYTICAL PROPERTIES OF THE POLARIZABILITY FUNCTION 

We start by defining an extension to complex values of q of the function x° defined in Eq. (|4]). Let's define x'^il^P) 
for p > and complex q ^ x + iy by the following formula: 



mp 
'2^ 



m {q~ 2p){q + 2p) J |9 " 



2p\ 



2p\ 



+ i 



arctan 



2p 



arctan 



y 



X — 2p 



(9) 



We first determine a maximal domain in the complex plane where the complex function is defined. For p > 



we define the complex domains Gi p, G2,p and G^^p to be respectively the set of all complex q 



iy such that 



JRq < —2p, —2p < ^q < 2p and SRg < 2p, respectively . One can check that the expression ^ defining 



converges to 



■mp 



as q — 0. We also consider L, the principal determination of the complex logarithm 

L{w) ~ log \w\ + i argw, 



where arg is the principal determination to the argument function, taking values between — tt and +7r. If q is in Gi^p 
we have the following representation: 



0/ \ '"-y 



mp m, [q — 2p) {q + 2p) 
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{L[- {q + 2p)]-L[- {q~2p)]}. 



We also have similar representations for g ^ in G2.p and in G^ p. It follows that x^{q.,p) is analytic in the domains 
Gp^i, excluding perhaps g = in Gp^2- Since x^{<liP) converges as g ^> 0, it follows that x^(q.,p) is analytic on the 
three domains Gp^i. 

Thus x^isiiP) is analytic in the whole complex plane, excluding the two branch cuts x = ±2p. One can also check 
that x^ can be continuously extended to the real points q — zL2p. It is also true that x'^iQ^p) can be continuously 
extended to any of the closures Gp^i, but those extensions don't match up at the boundaries (except at the real points 
q = ±2p) since x*^ has jump discontinuities there. 

We now consider x^ defined as in (jS]) using the complex x° instead. In order to prove that the integral defining 
exists and to study its analytical properties, (continuity or computation of complex integrals) it is necessary to use the 
classical theorems of Fubini and Lebesgue from integration theory. The key properties of the function g{p), (O, that 
will be needed are that g converges to zero exponentially as p — )■ oo and that g{p) — 0{p) as p — > 0. These properties 
will allow us to dominate products of g(p) and not too fast growing functions of p. In particular, if a function h{p) has 
polynomial growth as p — > oo and is continuous, or has a behavior like l/p or like logp as p — >■ , then the product 
g{p)h{p) will be Lebesgue-integrable and the necessary computations will make sense. 

We first consider for which points q = x -\- iy oi the complex plane does the integral in Eq. (jSj exist. First of all, 
since x"(0,p) — —mpliv^^ we have that x^(0) is well defined by Eq. If we fix a point go 7^ 0, formula ([9]) can 
be used and we can bound x''(90:?') by a- polynomial in p. In this way, we see that the function x^ is defined on the 
whole complex plane. 

Turning to the continuity of X"^j we can start with a point go 7^ and a sequence 7^ g„ — go- We consider the 
compact set K = {go} U {g„, n = 1, 2, . . .} so that ^ iiT. We can also use Eq. © to obtain a polynomial bound to 
|x(g,p)| uniformly on g G so that Lebesgue's dominated convergence theorem can be applied and the continuity 
of x^ at go 7^ can be established. 

To prove the continuity of x^ at g = 0, we bound x''('?j?') iii ^ slightly different way. We first consider q ~ x + iy 
with |g| < 1. If g £ G2,p \ {0}, we have 



x"{q,p)^-^ + £^{q-2p){q + 2p)i;p{q,) 



where 



V'p(<?) 



L{2p + g) - Li2p - g) 



L is the principal determination of the complex logarithm as before, and the function ipp can be extended to an 
analytic function in the whole domain G2,p. One then checks that 



V'p(g) = 2L ( 1 + ^ 
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Figure 1: The electron gas polarizability at zero temperature (dashed hne), for a value qi — 1, as a function of qr- All 
magnitudes are in units of the electron mass. For the chosen parameters (see the text), we have 2p — 0.00934 . Also shown for 
comparison is the finite-temperature polarizability (solid line), when T = O.OlTi?. 



We can check the above equality by expanding both functions in power series for small ^ and applying analytic 
continuation to G2,p- Then we have the estimate 

< |log|2p + g|| + |log2p| + |, 

which allow us to bound Ix^l'ZiP)! by the sum of a polynomial on p and a term which is the product of a polynomial 
on p and | log2p|. 

On the other hand, ii q E Gi.p U 6*2, p, the estimate of x'^l'ZiP) a-s in ^ gives the sum of two terms. The first term 
involves the logarithmic part of ([9]) and gives, after some inequalities, a polynomial bound. The second term involves 
the "arctan'' part and, using the fact that q £ Gi^p U G^^p gives a bound by a polynomial in p multiplied by 1/p. 

In any case, if \q\ < 1, x'^iQiP) is bounded by terms that are polynomial on p, logarithmic or 1/p, uniformly on q. 
Therefore, the product g(p)x'^{QiP) is dominated by a Lebesgue-integrable function and is defined and continuous 
in the whole complex plane. 

Thus, at finite temperature the above mentioned discontinuities disappear, and the function x^io) becomes con- 
tinuous everywhere. This is clearly illustrated in Fig. [U where we have plotted the real part of x°('?;P) (dashed 
line) . We have chosen p as the Fermi momentum of an electron gas with rj = 3 , where rg is the mean interelectronic 
distance in units of the Bohr radius, as usual. We also plot (solid line) the real part of x^il) for the same density, 
and a temperature T = O.OlTp (here, Tp is the Fermi temperature). As is apparent from this figure, the discontinuity 
of x^^il^P) around 5ft(q) = 2p disappears when the temperature is non zero. 

It turns out that x^ is non analytic at any point of the complex plane. To show this, we can compute the integral 



dqx^iq)^ dq dp g{p)x\q,p), (10) 

R JR. Jo 

where i? is a rectangle. If the rectangle is contained in the first quadrant, standard estimates on |x°('ZjP)| show that 
the double integral of x'^il^P) is well defined, the needed change of the order of integration is possible, and the usual 
procedures in complex analysis show that the value of (jlOp is not zero. The computation can be reproduced in all the 
other quadrants and we finally arrive to the conclusion that is not analytic at any point. 

We give the details of the computation in the first quadrant as an illustrative example. Suppose that the rectangle 
R is defined by the vertices q — xq + iyo and q = Xi + iyi. 



We introduce the notations 



X°(C2+) = lim^x"(2p-e + %,p), 



and the definition 



X°(C2-) = linixO(2p + e + %,p), (11) 



AX"(»,P) . x'iCi, - x"(C-,- ) = (12) 
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By deforming the contour integral, we can write 



Therefore 



which is different from zero. 



C Jyo 



rxi/2 ny-i 

dqx^{q)^-il dpg{p) dyAx°(y,p), 

R Jxo/2 Jyo 



IV. COMPUTATION OF THE POTENTIAL 



We now return to the computation of V{r), as given by Eq. 
First, let us calculate 



(13) 



(14) 



dqqcxp{iqr)x {q) = / dpg{p) / dqqexp{iqr)x {q,p), 
c Jo Jc 



(15) 



where C is the contour in figure [51 

In order to see that the above integral is well defined, to be able to interchange the order of integration and to 
compute the limit when the width of the contours in figure [2] tends to zero, it is necessary to estimate x^{q): and take 
into account that the imaginary part of q in the above contours is unbounded. 

We can keep the width of the contours in figure [5] equal to e < co- Then we can make the estimates for the full 
rectangles of width eg. It can be checked that the expression gexp(?gr)x°(q,p) admits a polynomial bound in p 
uniformly on q inside those rectangles if y < and that it can be bounded by the product of a polynomial on p by 
the exponential exp( — |g|r/2) \i y > \x\. With those estimates, it can be seen that, defining 

! d\q\ \qe^'''^9{p)x°{q,P)\. 
Jc 

it follows that p I(jp) < +oo, and the order of integration can be interchanged. Indeed, the expression 



-lip) 



dqqe^'^''x°{q,p) 



can be dominated with a Lebesgue-integrable function of p, so that the limit when e — > can be taken inside the 
integral. 

We now compute the value of (|T5|) . In the following, we leave aside the contribution from the Debye pole, which 
lies on the imaginary axis. At zero temperature, its position is given by qu = ±i(pjmA)^/^/7r. The pole contribution 
can be easily incorporated in our calculations via the residue theorem, giving rise to an exponentially damped term, 
which does not appreciably modify the results for the long range behavior we are interested here. Moreover, at finite 
temperature the pole is located at higher positions on the imaginary axis [loj , therefore giving an even smaller 
contribution. 

Under this approximation, we can write 



dqqexp{iqr)x°{q,p) 

L+ dqqexp{iqr)x°{q,p) - L- dqqexp{iqr)x°{q,p) 
L+ dqqexp{iqr)x°{q,p) - L- dqqexp{iqr)x° {q,p)dq 

Let us define, in analogy to Eq. ([TT|) 



(16) 
(17) 



X°{C%) ^ limx°(-2p- 



X^Cl2p) ^ limx^-2p 



-iqt,p) 
-iqi,p) 



(18) 
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y = Jmq 



-2p 



2p 



Figure 2: Analytically-continued structure of the zero-temperature polarizability. Represented schematically is the contour 
path used for integration. 



Due to Schwartz's principle, one has: 



which implies that 



Substitution on ([16]) gives, after some algebra 

/ dqqexp{iqr)x°iq,p) ^ ~2i dp exp{-yr)Re[qexp{i2pr) Ax"iy,p) 
Jc JO 

Therefore 

/ dqqexp{iqr)x^{q) = -2i dpg{p) dy exp{~yr) 
Jc Jo Jo 

Re [{2p + i) exp {i2pr) Ax°{y,p)] 
By repeating the above procedure, one can obtain the result: 



(19) 



(20) 



(21) 



(22) 



dq q exp{iqr) [x^ (q)] = -4« / dpg{p) dy exp{-yr) 

JO JO 

Re [{2p + iy) exp {i2pr) Ax'' {y , p)x^ {2p + iy)] 
which can be easily generalized to the following formula 



I dqqcxp{iqr) [x^{q)Y = -2ni dpg{p) dy exp{-yr) 

JC JO JO 



Re 



{2p + iy) exp {i2pr) Ax°{y,p) [x^{2p + iy)] 



(23) 



(24) 



Again, the necessary estimates on the above functions, in order to justify all the steps, can be done in a direct 
albeit long procedure. 

With the help of the previous equations, we can finally proceed with Eq. U). 



I{r) = Im dq i^^HM ^ dq ■ " '^^P^^'^^^ 



K?) 



Xx^iq) 



(25) 
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To this end, we make an expansion of the denominator in powers of A. Using the result of Eq. (1^^ we obtain, after 
a straightforward calculation: 



/(r) 

Here 



-2 / dpg{p) / dy exp{-yr)Re 
Jo JO 



AAx°(y,p) 

q exp (iZpr) 



u{qf 



(26) 



q = 2p + iy (27) 



We can now use the above result to obtain an approximate expression for I{r) as r — s- oo. In this case, due to the 
fast-decaying exponential, is is enough to consider only small values of y. To the leading order in y we have, then: 

Ax (y,p) ^ - — 

v{2p + iy) ~ v{2p) (28) 



and one easily obtains 



\m f°° pcos{2pr) 
7rr^ Jo i'(2p)^ 



m^^l dpgipf-^^^ (29) 



This formula is valid, for an arbitrary temperature, at sufficiently large distances and is specially suited, in contrast 
to the initial expression Eq. ([T|), to low temperatures. Indeed, as T — the function g{p) is strongly peaked around 
the Fermi momentum pp and allows for a fast convergence in the above expression. Within this limit, therefore, we 
can make further approximations, namely: 

/(r)^^-^/ dpg{p)cos{2pr) (30) 



7rr2 v{2pFy Jo 
After some algebra, one arrives to the final result 



rrr)-g^' ""'^ cos(2ppr) 



PF 



In this way, we have obtained an oscillatory behavior for V{r) at finite temperature. Whenever "^^"^J^ ^ 1, one 

can approximate sinh (^^^^^^ — \ exp ^ ^ i giving raise to an exponential suppression of Friedel oscillations at 

finite temperature. This is a very well-known result, which can be obtained using e.g. some modification of Lighthill's 
method to finite temperature (see, for example (ill. [T2|). Such methods, however, immediately raise the question of 
what the possible origin of the oscillations is, as discussed in the introduction. Our method is based on the integration 
of %''((?, p) on the complex plane, by interchanging the limits of integration (see p5l) V In the next section we give an 
intuitive description of this procedure, leading to pip . 



V. DISCUSSION 



In this paper, we have revisited the method to obtain the screened potential on an electron gas at finite temperature 
using the properties of the zero-temperature polarizability. We analyzed with detail the mathematical properties of 
this function on the complex plane. The formula we obtained for the potential in the low-temperature regime coincides 
with previous results in the literature, giving an oscillatory function, which is damped as temperature increases. The 
existence of oscillations at non-zero, but sufficiently low temperatures, shows that the limit T is smoothly reached 
in the model, and does not represent a pathology. However, the explanation based on the Kohn singularity fails to 
account for the persistence of the oscillations. In fact, the electron distribution is smeared out over a range of the 
order ~ T in energies, and therefore the Kohn singularity disappears. On the other hand, for T > the function X"^('?) 
becomes continuous, so that the argument we used at zero temperature does not apply to explain Friedel oscillations 
at finite temperature. 

In fact, as we have learned from the previous section, the function 'X^ [q) is non-analytical everywhere on the 
complex plane. From an intuitive point of view, we can regard Eq. ([3]) as a superposition of a family (as p varies) 
of functions x'^((7,p), each one having discontinuities at gr = ±2p. This sum is spread out, due to the weight function 
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g{p). over a range ^ T in energies. By interchanging the order of the integration, as in (jlSp . we obtain a superposition 
of oscillatory terms, Eq. (|5D)) . which results in (|3ip . 

In our calculations, we have considered only a qualitative discussion based on the RPA. It would be interesting to 
go beyond this approximation, by introducing a local field correction. Unfortunately, to our knowledge, there is no 
analytical formula which accounts for these effects in the case of finite temperature. 
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